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Abstract 

This article describes a normal form algorithm for the Brieskorn lattice of an iso- 
lated hypersurface singularity. It is the basis of efficient algorithms to compute the 
Bernstein-Sato polynomial, the complex monodromy, and Hodge-theoretic invari- 
ants of the singularity such as the spectral pairs and good bases of the Brieskorn 
lattice. The algorithm is a variant of Buchberger's normal form algorithm for power 
series rings using the idea of partial standard bases and adic convergence replacing 
termination. 
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1 Introduction 



Isolated hypersurface singularities form the simplest class of singularities. 
Their intensive study in the past has led to a variety of invariants. The Milnor 
number is one of the simplest, and can easily be computed using standard basis 
methods. A finer invariant is the monodromy of the singularity. E. Brieskorn [1] 
developed the theoretical background for computing the complex monodromy. 
He gave an ad hoc definition of an object H" , later called the Brieskorn lattice. 
Its great importance was a priori not clear. The complex monodromy can be 
expressed in terms of the differential structure of the Brieskorn lattice. 
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The finest known invariants come from a mixed Hodge structure associated 
to an isolated hypersurface singularity. The notion of a mixed Hodge struc- 
ture was introduced by P. Deligne [2] as a generalization of the classical Hodge 
structure on the cohomology of a compact Kahler manifold. J.H.M. Steenbrink 
[3] defined this mixed Hodge structure in terms of resolutions of singularities. 
A.N. Varchenko [4] and later J. Scherk and J.H.M. Steenbrink [5] described this 
mixed Hodge structure in terms of the differential structure of the Brieskorn 
lattice. The mixed Hodge numbers correspond to the spectral pairs and de- 
termine the complex monodromy. The spectral pairs have a semicontinuity 
property [6] with respect to unfoldings of the singularity. 

Based on properties of the mixed Hodge structure, M. Saito [7] constructed two 
endomorphisms Aq and Ai of the Milnor algebra. These two endomorphisms 
determine the differential structure of the Brieskorn lattice and, immediately, 
the above invariants. 

The Bernstein-Sato polynomial is associated to a general complex polynomial 

[8] or convergent power series [9]. T. Oaku [10] presented the first algorithm 
to compute it in the global case. A new method by M. Noro [11] is impres- 
sively faster. In the isolated singularity case, B. Malgrangc [12] described the 
Bernstein-Sato polynomial in terms of the Brieskorn lattice and its close re- 
lation to the complex monodromy. 

In [13-16], we have developed algorithmic methods to compute all of the above 
invariants of isolated hypersurface singularities. There is an implementation 
[17,18] of these algorithms in the computer algebra system Singular [19]. Our 
algorithm to compute the complex monodromy is much faster and can com- 
pute much more difficult examples than Brieskorn's algorithm. Our algorithm 
to compute the local Bernstein-Sato polynomial is based on B. Malgrange's 
description in terms of the Brieskorn lattice. It is much faster than M. Noro's 
algorithm since computations in rings of differential operators are replaced by 
computing their action on power series rings. However, it is restricted to the 
isolated singularity case while M. Noro's algorithm works in the general global 
case. All our algorithms require the computation of a basis representation in 
the Brieskorn lattice. In [13, Sec. 10.2] this is done by a sequence of divisions 
by the Jacobian ideal which is, in general, very hard to compute. The subject 
of this article is a normal form algorithm for the Brieskorn lattice replacing 
this sequence of full divisions by a sequence of partial divisions. This new 
method turns out to be much more efficient. 

In the first section, we recall the definition and the main properties of the 
Brieskorn lattice. We introduce the formal Brieskorn lattice and describe it as 
a cokernel of a formal family of differential operators which is finite over the 
base. 
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In the second section, we consider such a formal family of differential operators 
in general. We describe a normal form algorithm to compute a presentation 
of the cokernel which is a finitely generated module over the formal power 
series ring in the parameters of the family. This algorithm is a variant of 
B. Buchberger's [20,21] normal form algorithm. There are three major differ- 
ences compared to the classical algorithm: 

(1) The polynomial ring is replaced by a formal power series ring. Termina- 
tion of the algorithm is replaced by adic convergence. 

(2) The standard basis is replaced by a partial standard basis, a set of power 
series which specializes to a standard basis. 

(3) There is only a module structure with respect to the parameters of the 
family and the partial standard basis is not finite. 

Although the algorithm does not terminate in general, it serves to compute 
exact results by using appropriate degree bounds. The algorithms in [13,16] 
implicitly contain such degree bounds to compute the above invariants of iso- 
lated hypersurface singularities. There are also a priori degree bounds in [15], 
but they are useless in practice. Essentially the double number of variables plus 
the double Milnor number is a degree bound that satisfies all requirements. 

In the third section, we demonstrate the power of our algorithm. We apply 
the Singular implementation [18] to examples from [13,16] and [11] and list 
the timings. 

Families which are finite over the base occur in many situations in algebraic 
geometry and singularity theory. For example, A. Friihbis-Kriiger [22] has de- 
veloped algorithms to compute moduli spaces and adjacencies of singularities 
based on the idea of partial standard bases. One can expect more applications 
of our methods in the future. 

We shall denote row vectors by a lower bar, column vectors by an upper bar, 
row indices by lower indices, and column indices by upper indices. 
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2 The formal Brieskorn lattice 



Let / : U >C be a holomorphic function on an open neighbourhood e 

U C C" of the origin. We choose a system of complex coordinates x — 
Xi, . . . , x„ at G C" and denote hj d = di, . . . ,dn = d^^^ ■ ■ ■ , dx„ the cor- 
responding derivatives such that the commutator of di and Xj is [(9^,3;^] = dij. 
We consider / as a germ of a holomorphic function at G C", which means 
that U can be arbitrarily small. This is equivalent to considering the conver- 
gent power series / e C{^}. We assume that /(O) = 0, and that the origin 
is an isolated critical point of /. This means that G ?7 is the only point 
with (9i(/)(0) = • • • = (9„(/)(0) = for some U, or, more algebraically, that 
{x)"^ C {d{f)) C (x) for some m > 1. The complex dimension 

// = dimc(C{x}/(9(/))) <oo 

is called the Milnor number. By the finite dctcrminacy theorem [23, Thm. 9.1.4], 
one can choose, in this coordinate system x such that / G C[a;] is a 

polynomial. 

Wc denote by Q* = i7* „ q the complex of germs of holomorphic differential 
forms at G C". Its elements arc differential forms with coefficients in the 
convergent power series ring C{x}. The Brieskorn lattice [1] is defined by 

//" = QVd/A dQ"-' 

and becomes a C{t}~module by setting 

t-[u;]^ [M (1) 

for [a;] G H". By M. Sebastiani [24], H" is a free C{t}-module of rank /x. We 
denote by fl the /x-dimensional C-vector space 

Q = QVd/AQ"-^ ^C{x}/(9(/)). 

The operators d and d/ = d/ A ■ define two exact sequences. 

Lemma 1 (Poinccire Lemma) 

>C{x}^^n^^^ ^0 

is an exact sequence of C-vector spaces. 

Since completion is exact. Lemma 1 remains valid when replacing Q* by its 
(a;)-adic completion ^2*. The elements of Q* are differential forms with coeffi- 
cients in the formal power series ring C|x]. 
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Lemma 2 (De Rham Lemma) 

— >c{x}^n^^---^n- — >n — 

is an exact sequence of C{x} -modules. 

Also Lemma 2 remains valid when replacing Q' by Q*. From Lemma 1 and 2 
follows that one can define a C-linear operator s on H" by 

s-[drj]^[dfAri] (2) 

for [d?7] G H". From Lemma 1 follows that s is injective. The image of s is 
sH" = df A ^T/ df A df]"-2 and hence 

H"/sH" = n. 

Also s defines a module structure on H" over a power series ring. This power 
series ring is the ring 

of microlocal operators with constant coefficients and, by F. Pham [25], H" 
is a free C{{s}}-module of rank ji. From the definitions of t and s follows 
immediately that the commutator of t and s is 

We define a C-hnear operator dg on the localization i?"<8)c{{s}} C{{s}}[s~^] by 

t = s^a,. (3) 

Then i is a differential operator on H" with respect to the C{{s}}-structure. 
There is also C-hnear operator dt on the locahzation H" ®c{t}'^{'t}[t~^\ defined 

by 

s^d;\ (4) 

Then the commutator of dt and t is [dt, t\ — l and hence dt is a derivative by 
t. 

Definition 3 

(1) We call the topology induced by the (x) -adic topology on fl"' on the quo- 
tient H" the {x)-adic topology on H". 

(2) We call the completion H" of H" with respect to the {x)~adic topology 
the formal Brieskorn lattice. 
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Lemma 4 There is a natural isomorphism 

H" = Q^'/df A dQ''-\ 

Proof. By definition, 

H" = lim(Q"/((x)*^Q" + d/ A dQ"-2)). 

k 

Since e C" is a critical point of /, C (x) and hence 

d/ A d[{x)^Q''-^) C (x)'^fi", 
d/A d((x)'=0"-2) c 

Since / {x)^^r = tL"" / {x)''^r , tliis implies that 

^''/({xfQr + d/ A dQ"-2) = Qrl[{xfQr + d/ a dQ"-^) 

defines a natural isomorphism of inverse systems. Hence, 

H" = \^(pr / ({x)^QP + d/ A dQ"-2)) = Q"/ d/ A dOP-^. □ 



The following theorem [1, Prop. 3.3] is essential for Bricskorn's algorithm to 
compute the complex monodromy, which is based on the C{t}-structure of 
the Brieskorn lattice. 

Theorem 5 The {t)-adic and {x)-adic topology on H" coincide. In particu- 
lar, the (t) -adic completion of H" is naturally isomorphic to H" and H" is a 
free 'C\t\-module of rank /x. 

The C{{s}}-structure of the Brieskorn lattice is more algebraic and, therefore, 
more appropriate for computational purposes. The following proposition [16, 
Prop. 7] is the analogue of theorem 5 for the C{{s}}-structure, but it is much 
easier to prove. 

Proposition 6 The {s)-adic and {x)-adic topology on H" coincide. In par- 
ticular, the (s) -adic completion of H" is naturally isomorphic to H" and H" 
is a free C[s]-moo?uZe of rank /i. 

Proof. We denote 

dx — dxi A • • • A dxn, 

dx^ — dxi A • • • A dxi-i A dxi+i A • • • A dXn 
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ioT 1 < i < n. Let 

gdiif) dx] e dx+dfA dQ''-^) /df A dQ''-^ C H" 

for some A; > 1. By (2), 

[gdiif) dx] = [i-iy^'dfA^gd^)' 
(-1)^+M(^d^)] 
dx 

and hence, by induction, 

{{d{f)f''dx+ df A d^r-^)/ df A dQ"-2 c s^H". 

Since e C" is an isolated critical point of /, (x)"* C C {x) for some 

m > 1 and hence 



\n-2 



({xf'''^dx+ df A d^r-^)/df A dU'- 

C {{d{f)f^dx+ df A df)"-2)/d/A dJ^'^-'. 



This imphes that 



[{xf^'^dx^ df A dQ"-')/d/ A dOr-^ C 5*^/7" 

C ((x)'=dx+ d/A dQ"-')/d/A dQ 

Hence, the (s)-adic and (x)-adic topology on H" coincide. □ 



,n-2 



Note that the formal Brieskorn lattice is a (t, s)-module in the sense of D. Bar- 
let [26,27]. The following proposition [16, Prop. 8] describes the C|s]-module 
H" as a quotient of the power scries ring C|s, x\. It will lead to a normal form 
algorithm for H" in the next section. 

Proposition 7 dx induces a 'C\s\-isomorphism 

H" = Q"H/(d/ -sd)n--^lsl^Cls,x\/{d{f) - sd)Cls,xl 



Proof. Since 

df A dQJ'-'^ = ( d/ - s d) dQ"-^ C ( d/ - s d)f]"-^ Is] 
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and by Lemma 4 and (2), there is a natural C|s] -linear map 

Let u = Efc>o ^fes'^ e tl^'-'^lsj with ( d/ - s d)a; G 1]". Then d/ A Wfe+i = duok 
and hence, by (2), 

s[du;k+i] = [d/ Aa;fe+i] = [duk] G ^" 
for all A; > 0. In particular, [duo] e riA;>o s'^H" = {0} and hence, by Lemma 4, 

du;o e d/ A dQ"-2 = d( d/ A Q"-^). 
By Lemma 1, this implies that cuq G df2"~^ + d/ A and hence 

(d/-sd)a;= d/Aa;o e d/A dQ"-^. 

This shows that 

(d/ - sd)Q"-^|s] n = d/ A dQ"-2 

and hence, by Lemma 4, that is injective. By Lemma 1, dO**"^ = O", and 
hence is surjective. 

For 77 = Er=i(-1)'+Vid^e ^^"-'W, 

n 

( d/ - s d)77 = Y^{di{f)gi - sdi{gi)) dx = - sd)gdx. 

i=l 

Hence, dx induces a C|s] -isomorphism 

^^isl/(d/ -sd)n-^ls\^Cls,x\/{d{f) - sd)Cls,x\ . □ 



Proposition 7 is the starting point for more general considerations in the next 
section. 



3 Formal families of differential operators 



Let K be a computable field and F = Fi, . . . , G IK[s, x\{d) a formal family 
of differential operators where x = xi, . . . ,Xn, d = di, . . . ,0^ = dx^, ■ ■ ■ i dx„, 
and s — si,...,Sm- Note that the elements of K[s, are polynomial in d. 
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The brackets (•) indicate that the commutator [a;j,9j] = is not zero. We 
want to compute the cokernel H — K|s, x]/ x} of the K|s]-hnear map 

Kb, xf Kls, xj H > . 

We assume that the speciahzation 

/ = /i,...,/r = Z(s = o) eicy 

is independent of d and that C (/)K.|x| for some A; > 0. In particular, 

/X = dimK(Ky = dim^(H/{s)H) < oo 

and hence if is a finitely generated K|s]-module. Then there is a matrix 
D = {dP)j e such that 

F^f-sD. 



Our considerations are motivated by the following special case. 

Remark 8 By Proposition 7, for K — C, m — 1, r — n, f — d{f), and 

TJ r-^ TtII 

tl =C[sl -H 

is the formal Brieskorn lattice. 

Let <x be a local degree ordering with respect to a weighted degree deg^ on 
the set of monomials {x- 1 a e N"} = of ]K|x]. This means that 

W — 

is a semigroup homomorphism with degj.(xj) < 0, and that <^ is a semigroup 
ordering such that 

deg^ix-) < deg^{x^) ^ x- <xX^. 

The support of p = Y^aVaXr- £ K|x] is defined by supp(p) = {a € | 7^ 
0}. We denote the leading exponent, resp. leading term, with respect to <x 
by lexp, resp. lead. This means that 

lexp(p) = maxsupp(p), 

lead(p) = pie.p(p)a;''^^P(^) 

for p — PaXr- e and that 

lead(P) = {lead(p) |p e P} 
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for a subset P C IK[a;]. Note that the maximum exists by Dickson's Lemma 
[28, Lem. 1.2.6]. The weighted degree deg^. extends to by setting 

deg^(p) = deglexp(p) 

for p e IK|x]. Since [di,Xi] — 1, deg^. extends to by setting 

deg^(9i) = - deg^(a;i) > 0. 



Let g = gi, . . . ^gihe & standard basis of (/)]K|[x]. This means that ^ gi & 
(/)%] and 

lead((/)Ky) = (lead(£))Ky (5) 

which implies that 

imix] = {g)Klx] (6) 

by the division theorem. Let 

m = ("T.i)j=l,...,^ = (^-)x^^(lead{g))K|y| 

be increasingly ordered with respect to <^. Then 

= (m)K e (lead(£))K|x] (7) 

and hence, by (5) and (6), = {rn)K Q) (gf)K|x]. Then m represents a 

K-basis of 

H/{s)H = Klxj/igMx] = {m)K (8) 

and, by Nakayama's Lemma, m represents a minimal set of ]K|s]-generators 
of H. Note that if H is free then it is free of rank Let U — {ui)i G IK|x]''^' 
be a matrix such that 

g^fu. 

Remark 9 // / G K[x\ then one can compute g and U with coefficients in 
K.[x\ by hazard's method based on Buchberger's standard basis algorithm [29, 
Lem. 1.7] and homogenization. In general, the power series fi G K[a;] can 

he represented by generating functions s-K and Buchberger's standard 

basis algorithm with respect to a local degree ordering computes such generating 
functions for g and U. 

Let <s be a local degree ordering with respect to a weighted degree deg^ on 
the monomials of K|s]. Let 

<= i<s, <x) 

be the block ordering of <s and <x on the monomials of x] and 

deg = deg^ + deg^ 
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the sum of the weighted degrees deg^ and deg^.. This means that 

s^'x^ < s^'x^' ^ <, s^' V [s^ = s^' A x^ x^") 

and 

deg(s%^) = deg,(s^) + deg^(x^). 

As before, we denote the leading exponent, resp. leading term, with respect 
to < by lexp, resp. lead, and extend deg to ]K|s, Note that < is not a 

degree ordering with respect to deg. This means that 

deg lead ^ deg . 

We denote by the leading exponent, resp. leading term, with respect to the 
partial ordering <s on K|s, x] by lexp^, resp. leads, and the partial degree <s 
on x| by deg^. This means that 

lexp^(p) = max7rN'"(supp(p)), 

lead,(p)=pie.p^(p)si^^P^(^), 
degs(p) = deglexp^(p) 

for p = J2aPaS- G IK|s, where ir^m : N™ x N*^ is the canonical pro- 
jection. We denote by mindeg, resp. max deg, the minimal, resp. the maximal, 
degree of the components of a vector or a matrix. 

Let 

G = Gi,...,Gi = FU = g- sDU. 

In the special fibre s = 0, G induces the standard basis g. We call G_ a partial 
standard basis of the formal family F_. 

The following example is taken from [16, Sec. 8]. 

Example 10 Let K = C, F = d{f) — sd as in Remark 8, and f = x^ + x'^y'^ + 
y^. Note that f defines a T2^5,5 singularity at the origin. Let <(x,y) be the local 
degree ordering with deg(x) = deg(y) = —1 and x > y. Then one computes 

I = 2xy^ + 5x^ 2x'^y + 5^/^ 

g = 2x^y + 57/^ 2xy'^ + 5x^ bx^ - 5^/^ lOy^ + 25x^7/^ 
/X = 11, 

m = ^/^ y^, l/^ |/^ xy, y, x^ a;^ x^, X, 1, 



u 



1 X —2xy 
1 -y 2y2 + 
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and hence 

G = 2x^1) + — sdy, 2xy^ + bx'^ — sdx, 

5x^ - 5y^ - sxd^ + sydy, lOy*^ + 25x^y'^ + 2sxydx - s{2y'^ + bx^)dy. 

We denote by Fx^" = {FiX-)i^a the generators of the C|[s]-module (£)K|[s, xj . 

Lemma 11 H is a free Kls}-module if and only if Gx^" is a standard basis 
of the -module (Z)K|s, x] . 

Proof. By (8) and Nakayama's Lemma, m represents a minimal set of gen- 
erators of H. Since H = K|s,x]/(F)]K|s,x], 

= {m)Kis} + 

and H is free if and only if 

n (m)Kls] = 0. 

By (5) and (7), this is equivalent to 

lead((F)Kls,a;]) = (lead(^))Kls, = (lead(Ga;^"))Kls]. □ 

By Proposition 7, 

Cls,x]/{d{f)-sd)Cls,x]^cM H" 

is a free C|s]-module of rank /i. We shall now give an elementary proof of this 
fact. 

Proposition 12 IfK — C and F_ — d{f) — sd as in Remark 8 then H is a 
free ClsJ-module of rank /i. 

Proof. Let 7^ p e (zil)C|sl n - sa)C|s,x]. Then lead(p) e (m>C|s] 
and p = {d{f) + sd)q for some q G C[s, x] with maximal maxdeg3(g). By (5) 
and (7), this implies that 9(/)leads(g) = and hence, by Lemma 2, we may 
assume that there are l<i<j<n, k>0, and r e C|x] such that 

lead,(^) = sV(9,(/)e,- - 9,(/)e,). 

This implies that 

91ead,(g) = s''d{f)(dj{r)ei - 9,(r)e,) 
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and hence 



This is a contradiction to the maximahty of meixdegg{q). Hence, 



(m)Cls] n {d{f) - sd)Cls,x] = 



and H is free. □ 



Our aim is now to define a filtration V = {Vk)k<o on K|s, xj by K|s]~modules 
which is 

(1) a basis of the {s,x)-a,dic topology on ]K|s, ^J, 

(2) compatible with reduction with respect to the partial standard basis G_, 

(3) mapped by tth onto the basis {{s)^ H)k>q of the (s)-adic topology on H. 

This will lead to a normal form algorithm for H^. 

For a given weighted degree deg^,, let the weighted degree deg^ be such that 



Let V = {Vk)k<o be the strictly increasing filtration on x] by IK|s]- 
modules 



deg(sj) < mindeg(7Ti) + mindeg(x) — maxdeg(d-'). 



(9) 



Let the strictly increasing sequence N — {Nk)k<o be defined by 



Nk — — -fC mindeg(s) 



mindeg(x) + maxdeg(D). 



(10) 




(11) 



Remcirk 13 For F = d{f) — sd as in Remark 8, we can choose 



deg(s) = min deg(m) + 2 min deg(a;), 
Nk — —K deg(s) — 2 mindeg(x). 



Example 14 In example 10, deg(s) = —7 and Nk — IK + 2. 



The following proposition is a generalization of [16, Lem. 10]. 



Proposition 15 



(1) V = (Vft:)_ft:<o is a basis of the {s,x)-adic topology. 

(2) Iflead{s^GkX^) G Vk then also s^GkX^ E Vk- 

(3) MVk) = {sy^'H. 
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Proof. 

(1) This follows from (10) and (11). 

(2) Since g is a, standard basis, 

mindeg(m) +mindeg(x) < mindeg(gf) (12) 
and hence, by (9), 

deg{sDuk) < max{deg{sj(Fuk) | 1 < j < rri} 

< max{deg(sj) + maxdeg{(f) + maxdeg(Mfc) | 1 < j < m} 

< max{deg(sj) + maxdcg(c|'') | 1 < j < m} 

< min deg(r?i) + min deg(a;) 

< mindeg(£) < deg{gk) 
= deglead(g'fc). 

Since Sr^kX- = s-{gk — sDuk)x-, this implies that 

deglead(s^G'fcX^) = deg(s^lead(yfe)x^) = deg(s^G'feX^). 

Hence, the claim follows from (11). 

(3) Let Q ^ p E Vk and s-p^ = lead^(]?) with maximal |a| < —K for fixed p 
mod (Z)K|s,x]. Then, by (9), 

deg(pa) = deg(s-j9^) - deg(s^) 

< —{K + |a|) mindeg(s) — mindeg(x) + maxdeg(£') 

< min deg(s) — min deg(a;) + maxdeg(L') 

< mindeg(m) 

and hence, by (7), Pa G ((7)K|a;]. By the division theorem, there is a 
q e with Pa = gq and lea.d{pa) > lea.d{gjq^) for all j and hence 

maxdeg(g) < deg(pa) - mindeg(£). (13) 

Then 

Pa^9q^fUq = sDUq mod (F)K|s,xl (14) 
and hence, by (9), (12), and (13) 

deg{sDUq) < max deg (sZ^) + maxdeg(?7) + maxdeg(g) 

< maxdeg(sD) + maxdeg(g) 

< maxdeg(sD) — mindeg(5') + deg{pa) 

< maxdeg(sD) — mindeg(2Ti) — mindeg(x) + deg(pa) 

< max{deg(sj) + maxdeg(d-') | 1 < j < m} 
— min deg(7Ti) — min deg(x) + deg(pa) 

< deg(pa:)- 
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Hence, by (14), 



p' ^p- lead,(p) + s^sDUq = p mod (Z)K|s, xj 

with deg(p') < deg(p) < and lead5(p') <s leaded?). This contradicts 
to the maximahty of \a\ and hence p e {s)~^ + (£)IK|s, □ 



Proposition 15 leads to the following normal form algorithm. 
Algorithm 1 

proc NF(p e K|s, x], < 0) 
i/p e then q :— p 

else i/deglead(p) < Nk or lead5(p) e {s)~^ then q :— leads(p) 
else q := 
r := p — q 

ifr — then return r e K|s,a;], a e K|s,^]', q e IK[s,^] 
«/lead(r) e (lead(£)) then 

j min{i | lead(r) G (lead(c/i))} 

r,a,q' := NF(r - i^^. - i^) 

— — I lead(r) — 

else 

r',a, q' := NF(r - lead(r), /sT) 
r := lead(r) + r' 
q:^q + q' 

return r e K|s, x] , a e K|s, x]', e x] . 

The input of the algorithm NF is a power series p G K|s, and an integer K < 
0, the output is a power series r G ]K|s, x], a column vector a with coefficients 
in ]K|s, x], and a power series q G ]K|s, x]. We denote the components of NF 

by 

(NFi(p, K), NF,{p, K), NFsip, K)) = (r, a, g) = NF(p, X) 
for p G K[s, x] and i^T < 0. 

Example 16 In example 10 using 14, one computes 

NFi(/m, -2) = rn{Ao + sAi) 
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where Ao,Ai e C^^"^^^ such that 
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Figure 1 illustrates a reduction step in NF. The 1-dimensional K- vector space 
spanned by a monomial in (m)]K[s], resp. in (lead((7))K|s, x], is depicted by a 
big, resp. small, bullet. The monomial at the tail of the arrow is replaced by a 
power series with support above the dotted line meeting the head of the arrow. 
The K|s]-submodule Vk generated by the monomials above the dashed line 
is invariant with respect to such a reduction step. 



S A 



K 



Fig. 1. A reduction step in NF 



K 



X 



m)Ky : (lead(£))K[s,2;] 



Lemma 17 NF terminates. 
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Proof. For fixed leading exponent lexp^(p) with respect to s, the leading 
term lead(7;) is strictly decreasing with weighted degree deglead(p) > N^- 
Since there are only finitely many monomials with fixed lexp^ and deg > Nx, 
lexp^(p) decreases after finitely many steps. Since <= (<5, <^) is a block 
ordering and <s is a degree ordering, this implies that p e (s) ~^ after finitely 
many steps. □ 

The following lemma states that NFi(-,L) is a reduced normal form modulo 
Vl with NFi(y^, L) C for L < X < 0. 

Lemma 18 Let L < K < Q, p e Vk, and {r,a, q) = NF(p, L). Then 

(1) p^Ga + r + q, 

(2) a* e IK[s, x] with lead(G'ia*) < lead(p) for i = 1,. . . ,1, 

(3) r E (m) ®L<\a\<K^s^ with lead(r) < lead(p), 

(4) qeVL. 

(5) Ifp e (m) 0L<|a|<i^ Ks^ then (r, a, q) = (p, 0, 0). 

(6) Ifp = p' mod Vl then r — r'. 

Proof. By Proposition 15.2, NF preserves the condition p e Vk- Hence, the 
claim follows immediately from the definition of NF. □ 

By Proposition 15.1, V is a basis of the (s, x)-adic topology and, by Lemma 
18, NFi(y^,L) C Vk and m^iVK^L) C Vl for L < K < 0. Since Kls,xj is 
complete with respect to the (s, a;)-adic topology. 



and hence NF induces a reduced normal form on K.|s, x] as follows. 
Definition 19 Let K = (ii'j)j>o he a strictly decreasing sequence and 



for p e ]K|s, x\ where Po—p and ri, ai,Pi+i — NF(pj, Kj) for i >0. 
Note that NF depends on the choice of the sequence K. 
Lemma 20 Let (r,a) = NF(p). Then 

(1) p = Ga + r, 

(2) e x] with lead(Gja^) < lead(p) for i = 1, . . . ,1, 



Kls,x]^lim{Kls,x]/VK) 



K 
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(3) r e (m)IK|s] with lcad(r) < lca.d{p). 

(4) Ifpe (m)K|s| then {r,a) = (p,0). 



Proof. This follows immediately from Proposition 15.1 and Lemma 18. □ 



The following proposition describes NFi as a map of C[s]-modules. 
Proposition 21 NFi is a K|s]-/mear map 

Hs, ^ — ^ {m)Kls} 




H 



with NF^ = NFi. In particular, if H is a free K|s]-moc?u/e then NFi induces 
the KlsJ-section 

NFi 

of the canonical projection hh with image NFi(if) = (m)IK[s]. This means 
that NFi is the m-basis representation. 

Proof. By definition of NF, NFi is K|s]-linear. By Lemma 20, NFi is a map 
over H. By Lemma 20.4, NFi is the identity on its image (i2i)]K|s]. □ 

If H is not a free K|s]-module then its relations can be computed as follows. 
By Proposition 21, 

H (m)Ky/((m)Ky n (Z)K|s,x]) 

= (m)ICy/NFi((F)Kb,xl) 
= (m)Ky/(NFi(Fa;^"))Ky 

and, in particular, by Lemma 18.3, 

H/isy^'H ^Kbi (m)Ky/((NFi(zx^"\y;,,x))Ky + (m)(s)-^Ky) 

where Fx^"\Vk is a finite subset of ]K|s, x]. 

Finally, we return to our starting point. Let K = C and F = d{f) — sd as in 
Remark 8. Then, by Proposition 7, H =c[s] H" is the formal Bricskorn lattice 
and, by Proposition 12, if is a free C|s]-module of rank /j,. We define the 
matrix e Clsj^''^' by 

mA^ = tm. 
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Then, by (3), A— + s'^d^ is the 221-basis representation of t. This means that 
there is a commutative diagram 



H"- 



By (1) and propsosition 21, 

mA^= NFi(/m). 

Example 22 In example 10 using 16, 

A^ = Ao + sA^ mod s^Cfsf ^^^^ 

By [16, Sec. 8], the non-diagonal terms of Ai and the terms of A— in s^C[sp^^"^"^ 
can be eliminated by transforming m to a good Clsj-basis of H" . Then 

A^^Ao + sA^ 

and Aq and Ai represent M. Saito's endomorphisms [7]. 

All the algorithms in [13-16] require the computation of the matrix A— for a 
C|s]-basis m of H". 



4 Examples and timings 



Algorithm 1 for the case of Remark 8 and the algorithms in [15,16,30] are 
implemented in the SINGULAR [19] library gmssing.lib [18]. Wc use this 
implementation on a Pentium III M 1 GHz machine with 512 MB of memory 
plus 1 GB of swap memory. For several polynomials / e C[x\ with isolated 
critical point at the origin, we compute 

(1) the local Bernstein-Sato polynomial, 

(2) the spectral pairs, and 

(3) M. Saito's endomorphisms and Ai. 

For i = 1, 2, 3, we denote by ti the corresponding computation time in seconds 
and by Ki the maximal K occurring in NF during the computation. By tx 
we denote the time in seconds needed to compute A— mod s^C|s]'*^'^. All 
computation times are rounded off. 

The local Bernstein-Sato polynomial at the origin for the examples in [11, 
Tab. 1] can be computed, each in less than one second. Table 1 shows the 
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timings for the examples in [11, Tab. 2]. By [15], it suffices to compute A— 
mod s'^"C|s]^^'^ where Ko = 2{fi + n— 1) in order to compute all of the above 
invariants. Table 2 shows the results for the examples in [13, Tab. 2] and that 
this a priori bound is useless in practice. 

Table 1 



Local Bernstein-Sato polynomial 6 for / = + y"^ _|_ ^"s + ^miym2^ms 



n 


6,6,6 


7,7,7 


7,7,7 


9,9,9 


6,6,7 


6,6,7 


6,6,7 


6,7,7 


m 


4,4,4 


2,2,2 


2,2,3 


3,3,3 


2,2,2 


3,3,3 


4,4,4 


2,2,2 




125 


167 


216 


512 


138 


150 


150 


152 


deg(6) 


13 


28 


17 


23 


60 


48 


52 


62 




3 


4 


2 


2 


4 


3 


3 


4 


h 


1 


5 


2 


50 


6 


1 


2 


6 



Table 2 

Spectral pairs and M. Saito's endomorphisms for [13, Tab. 2] 





f 




K2 




t2 






H 








x^y + x'^y^ + / 


16 


2 








4 





10 


34 


180 




x^ + x'^y^ + y'^ 


16 


2 








4 





2 


34 


39 


<1 


x^ + 2x^y^ + xy^ + y^ 


16 


2 





1 


4 


1 


54 


34 


495 


Q2,l 


x"^ + yz^ + x^y^ + y'^ 


15 


2 





1 


4 





4 


34 


58 


Q2,2 


x^ + yz^ + x^y^ + 


16 


2 








4 





1 


36 


20 


Si,i 


x^z + yz^ + x^y^ + y^ 


15 


2 








4 





1 


34 


43 


Sl,2 


x^z + yz^ + x^y^ + y''' 


16 


2 








4 








36 


27 


o# 
-^1,1 


+ yz^ + y^z + xy'^ 


15 


2 





1 


4 





18 


34 


687 


Q# 

'-'1,2 


x^z + yz^ + y^z + x^y^ 


16 


2 





1 


4 





2 


36 


252 


Ui,i 


x'^ + xz^ + xy^ + y^z^ 


15 


2 





1 


4 





2 


34 


66 




x^ + xz^ + xy^ + y'^z 


16 


2 








4 








36 


12 




x^y + z^y^ + z^ + y^ 


16 


2 








6 





2 


36 


12 


V* 

^1,1 


x^y + y^ + xz^ + yz^ 


16 


2 





1 


6 





1 


36 


10 



By Remark 8, (6) and (8), the coefficient of in the m-basis representation 
in E." is defined by a division by the ideal (9(/))C|£] where the output of 
the division for defines the input for the division for Therefore, 
the complexity of the data and the computation time increases rapidly with 
K. In [13, Sec. 10.2], we compute such a division by a sequence of weak 
normal form computations. Table 3 shows the time needed to compute A— 
mod s^C|sp^'^ for example 10 and increasing K using this method and the 
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Singular command division. The computation fails in degree K = S after 
more than one hour due to lack of memory. In the algorithm NF, the above 
sequence of full divisions is replaced by a sequence of partial divisions. Table 
4 shows the time needed to compute the same result using the algorithm NF. 
The situation is similar for other examples. 

Table 3 

A— for f = + x^y^ + using division 



K 


4 


5 


6 


7 


8 


iK 





1 


10 


283 


oo 



Table 4 

^— for f = + x'^y'^ + using NF 



K 


20 


40 


60 


80 


100 


120 


140 


160 


180 


200 


tK 


1 


2 


8 


18 


35 


64 


106 


165 


247 


354 
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